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Abstract 

The origin of galactic cosmic rays is one of the most interesting unsolved problems in astroparticle physics. Experi- 
mentally, the problem is attacked by a multi-disciplinary effort, namely by direct measurements of cosmic rays above 
the atmosphere, by air shower observations, and by the detection of TeV 7 rays. Recent experimental results are pre- 
sented and their implications on the contemporary understanding of the origin of galactic cosmic rays are discussed. 
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1. Introduction 

The Earth is permanently exposed to a vast flux of 
highly energetic particles, fully ionized atomic nuclei 
from outer space. The extraterrestrial origin of these 
particles has been demonstrated by V. Hess in 1912 
P] and he named the particles " Hohcnstrahlung" 
or "Ultrastrahlung" . In 1925 R. Millikan coined the 
term "Cosmic Rays". They have a threefold ori- 
gin. Particles with energies below 100 MeV originate 
from the Sun. Cosmic rays in narrower sense are par- 
ticles with energies from the 100 MeV domain up to 
energies beyond 10 20 eV. Up to several 10 GeV the 
flux of the particles observed is modulated by the 11- 
year cycle of the heliospheric magnetic fields. Parti- 
cles with energies below 10 17 to 10 18 eV are usually 
considered to be of galactic origin. A proton with an 
energy of 10 18 eV has a Larmor radius tl = 360 pc 
in the galactic magnetic field (B w 3 [iG). This ra- 
dius is comparable to the thickness of the galactic 
disc and illustrates that particles at the highest en- 
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ergies can not be magnetically bound to the Galaxy. 
Hence, they are considered of extragalactic origin. 

The energy density can be inferred from the mea- 
sured differential energy spectrum dN/ dE [2] 



4tt f E dN 
PE = -J JdE dE > 



(1) 



where 0c is the velocity of particles with energy 
E. For galactic cosmic rays the major contribution 
to the total energy density originates from parti- 
cles with energies around 1 GeV. Such particles 
are strongly influenced by the heliospheric mag- 
netic fields. Outside the heliosphere, in the local 
interstellar environment an energy density pj^ IS = 
1.1 eV/cm 3 is obtained (Ref. [3] p. 12). The param- 
eterization of the measured galactic cosmic-ray flux 
according to the poly-gonato model [4] results in a 
density p|? = 0.43 eV/cm 3 . This is the measurable 
energy density at Earth (for an average modulation 
parameter M = 750 MeV, see (1) in Ref. [4]). This 
implies that less than half of the energy flux can be 
registered directly at Earth. 

In this article we will give an overview on recent 
experimental results and their implications on the 
contemporary understanding of the origin of galactic 
cosmic rays. As space is limited here, the reader may 
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also consider further recent reviews by the author 
|4|5|6|7|8|. 

Progress in the understanding of the origin of 
galactic cosmic rays emerged mainly from observa- 
tions in three complementary disciplines. The direct 
measurement of cosmic-ray particles above the at- 
mosphere in outer space and on stratospheric bal- 
loons (Sect. [2]). At energies exceeding 10 15 eV the 
steeply falling cosmic-ray energy spectrum requires 
experiments with large detection areas exposed for 
long times, at present, only realized in ground based 
installations (Sect. QJ. With such detectors extensive 
air showers are detected, which originate in inter- 
actions of high-energy particles in the atmosphere 
(Sect. [3]). And, finally, the observation of TeV 7-rays 
(Sect.EJ). 



2. Direct measurements 

Detectors above the atmosphere allow a precise 
determination of the composition of cosmic rays and 
the measurement of energy spectra for individual el- 
ements, presently, up to about 10 14 eV. These find- 
ings are very important to understand the propaga- 
tion of the particles in the Galaxy. 

In the 100 MeV energy domain detailed studies 
of the isotopic composition of cosmic rays and the 
comparison to the local galactic abundance indicates 
that cosmic rays are accelerated out of a sample of 
well mixed interstellar matter [9] . The propagation 
time of the particles in the Galaxy has been deter- 
mined tober esc = (15.6 ± 1.6) -10 6 yr [10111112] . At 
GeV energies the elemental composition has been 
studied in detail, finding all elements of the periodic 
table in cosmic rays, see e.g. [6]. Some elements, as 
e.g. boron are produced during the propagation in 
the Galaxy through spallation of heavier elements 
such as carbon. Investigations of the abundance ra- 
tios of different elements (like the boron-to-carbon 
ratio) yield the column density of matter traversed 
by the particles in the Galaxy |10|13j . 

At present, several groups try to extend the direct 
measurements to higher energies towards the knee 
(at about 4.5 • 10 15 eV) in the all-particle spectrum, 
benefiting from long-duration balloon flights [15] . 
Among these are the ATIC p], CREAM [17], and 
TRACER [18j experiments. The latter is presently 
with an aperture of 5 m 2 sr the largest cosmic-ray 
detector for direct measurements. It comprises sev- 
eral electromagnetic detectors, namely an acrylic 
Cerenkov counter, a proportional tube array, and a 
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Fig. 1. Energy spectra for individual elements as measured 
by the TRACER, HEAO, and CRN experiments [H]. 



transition radiation detector, thus covering a large 
energy range from 10 GeV/n up to 10 5 GeV/n with 
single-element resolution. Spectra for elements from 
oxygen to iron as measured by TRACER during 
an Antarctic circumpolar flight in 2003 are com- 
piled in Fig.Q] [TJ] . The data extend previous results 
from the space experiments HEAO (satellite) and 
CRN (on the Space Shuttle) to higher energies. Cur- 
rently, the data of a flight from Sweden to Alaska in 
2006 are analyzed. TRACER has been flown with an 
upgraded electronics being able to record elements 
from boron to iron. The new data are expected to 
extend measurements of the boron-to-carbon ratio 
up to 10 5 GeV/n, thus, providing, unprecedented in- 
formation of the propagation properties approach- 
ing the knee in the energy spectrum. 

Frequently, the question arises whether the energy 
spectra of protons and helium have the same spec- 
tral index. Figure [2] shows a compilation of data, for 
details and references see Ref . [6] . From the figure it 
seems to be clear that the helium spectrum is indeed 
flatter as compared to protons. Due to spallation of 
nuclei during their propagation and the dependence 
of the interaction cross section oc A 2 / 2, one would ex- 
pect a slightly flatter spectrum for helium nuclei as 
compared to protons. Following [19] the difference 
should be of the order of A7 0.02. However, fits 
to the experimental data (lines in Fig. [2]) yield 7 P = 
-2.71 ± 0.02 and ~f H e = -2.64 ± 0.02 g], resulting 
in a difference A7 = 0.07. 
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Fig. 2. Energy spectra for protons (black) and helium nuclei 
(grey). For a list of experiments and references see Ref. [6]. 
The lines represent fits according to the poly-gonato model 

3. Extensive air showers 

When high-energy cosmic-ray particles penetrate 
the Earths atmosphere they interact and generate 
a cascade of secondary particles, the extensive air 
showers. Two types of experiments may be distin- 
guished to register air showers: installations measur- 
ing the longitudinal development of showers (or the 
depth of the shower maximum) in the atmosphere 
by observations of Cerenkov or fluorescence light 
and apparatus measuring the density (and energy) 
of secondary particles (electrons, muons, hadrons) 
at ground level. 

The shower energy is proportional to the total 
light collected or to the total number of particles 
recorded at observation level. More challenging is 
an estimation of the mass of the primary particle. 
It is either derived by a measurement of the depth 
of the shower maximum X max and the fact that the 
depth of the shower maximum for a primary particle 
with mass A relates to the depth of the maximum 
for proton induced showers as 

X max = X max ~ Xq In A, (2) 

where Xq = 36.7 g/cm 2 is the radiation length in air 
[20121) . Or, measuring the electron-to-muon ratio in 
showers. A Heitler model of hadronic showers [3T] 
yields the relation 

lg (iV e /i\g = <7- 0.065 In A (3) 

This illustrates the sensitivity of air shower exper- 
iments to In A. To measure the composition with a 
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Fig. 3. Number of electrons vs. number of muons at shower 
maximum for fully simulated showers (symbols). The lines 
represent predictions of a Heitler model: solid - constant 
mass for protons and iron nuclei J4), dashed — constant en- 
ergy JSJ, an d dotted - constant depth of the shower maxi- 
mum X m ax ©. 

resolution of 1 unit in In A the shower maximum has 
to be measured to an accuracy of about 37 g/cm 2 
(see or the Ne/N^ ratio has to be determined 
with an relative error around 16% (see ([3])). Due to 
the large intrinsic fluctuations in air showers, with 
existing experiments at most groups of elements can 
be reconstructed with A In A « 0.8 — 1. 

The detection principle is illustrated in Fig. [31 de- 
picting the number of electrons as function of the 
number of muons at shower maximum. The symbols 
represent fully simulated showers with discrete ener- 
gies in steps of half a decade, for details see Ref. [21] . 
The lines represent predictions of a Heitler model of 
hadronic air showers |21j . The solid lines are lines of 
constant mass 

N e \ A = 7.24 • A~°- 16 Nl- 16 (4) 

for primary protons and iron nuclei. The dashed 
lines correspond to a constant energy, following 

N e \ Eo = 5.30 ■ 10 7 (Eo/PeV) 1 - 37 A^ - 46 . (5) 

The sets of lines define a parallelogram giving the 
axes for energy and mass in the Ne-N^ plane as 
indicated by the arrows. The dotted lines represent 
a constant X max , separated by 100 g/cm 2 according 
to 
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Fig. 4. Number of electrons vs. number of muons for showers 
measured with the KASCADE experiment, the third dimen- 
sion is the muon production hight |22| . 

A similar plot is presented in Fig.Hl showing the 
N e -N^ plane for showers measured with the KAS- 
CADE experiment. The third dimension gives the 
production hight of muons [55] . In the N e -N^ plane 
light primary elements are expected in the upper left 
part of the figure. This is underlined by the measure- 
ments yielding in this area deeply penetrating show- 
ers. Attention should be payed when the two figures 
are compared: X max for the electromagnetic com- 
ponent (Fig. [3]) is compared to X max for the muonic 
component (E^ > 0.8 GeV, Fig.[4|. Fig.[3]represents 
N e and at the shower maximum, Fig.|4]gives val- 
ues measured at ground. For muons the lateral dis- 
tributions are integrated in the range interval from 
40 to 200 m only in Fig.g] 

By detecting two quantities (e.g. N e and or N e 
and X max ) , energy and mass of the primary particles 
can be deduced. The arrival direction of cosmic rays 
is determined by measurements of the arrival times 
of secondary particles at the detectors. 

4. Indirect measurements 

Air shower experiments extend the findings of di- 
rect measurements to highest energies. They provide 
information on the all-particle energy spectrum, the 
(average) mass composition, and, recently, informa- 
tion on energy spectra for elemental groups. These 
results contribute to the understanding of the origin 
of structures in the all-particle spectrum, like the 
knee or the second knee. 

The all-particle energy spectrum up to about 
10 18 eV is reasonably well known. The spectra ob- 
tained by different experiments agree well in shape 
and absolute flux taking a ±10% uncertainty in the 
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Fig. 5. Energy spectra for primary protons and helium nuclei, 
unfolded from air showers registered by KASCADE within 
three zenith angle intervals I24I . 



energy determination of the individual experiments 
into account [4] . The spectrum can be described by 
a power law dN/dE cx E 1 , with an spectral index 
changing at the knee (E^ = 4.5 PeV) from 7 = —2.7 
to 7 = -3.1. 

For the mean (logarithmic) mass of cosmic rays 
there seems to be a dependence on the technique 
applied for the air shower observations 23J. Two 
classes of experiments may be distinguished: detec- 
tors measuring the depth of the shower maximum 
and experiments measuring (the lateral distribution 
of) particles at ground level. This indicates incon- 
sistencies in the interaction models used to interpret 
the data. For example, introducing lower inelastic 
cross sections in the interaction model QGSJET 01 
and slightly increasing the elasticity, good agree- 
ment between the two classes of experiments can be 
achieved [23j . Interpreting the data with this model 
yields a consistent picture with a mean logarithmic 
mass smoothly continued from direct measurements 
up to about 10 15 eV, then exhibiting an increase up 
to energies around 10 17 eV. 

A big step forward in understanding the origin of 
the knee in the all-particle energy spectrum are the 
results of the KASCADE experiment. The latter 
comprises a 200 x 200 m 2 scintillator array to reg- 
ister the electromagnetic and muonic component 
[25], a hadron calorimeter equipped with liquid 
ionization chambers [26] . and a limited-streamer 
tube muon tracking detector [27j . The simultane- 
ous detection of different shower components allows 
detailed studies of hadronic interaction models 
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Fig. 6. Rayleigh amplitudes as function of energy as observed 
by different experiments |39| . 

[28I29I30I31I32I33I34] . The most consistent models 
are then applied to unfold energy spectra for five 
groups of elements from the measured air shower 
data [35]. The energy spectra of the light elements 
(protons and helium) exhibit a depression of the 
flux at high energies while the spectra of the heavier 
groups continue to follow power laws to higher ener- 
gies. The results obtained using different hadronic 
interaction models to interpret the data agree qual- 
itatively. However, there are differences in the ab- 
solute flux obtained, indicating inconsistencies in 
the interaction models. These results have recently 
been confirmed by an analysis taking also inclined 
showers into account [24] ■ As example, the spectra 
obtained for primary protons and helium nuclei for 
three different intervals of the cosmic-ray arrival di- 
rection are shown in Fig. [5] The different data sets 
are consistent with each other. 

Similar energy spectra for groups of elements have 
been obtained by the EAS-TOP [36], GRAPES-3 
[3"7] , and Tibet [35] experiments [5] . Overall, a con- 
sistent 'standard picture' evolves |4l6j : the energy 
spectra as obtained by direct measurements can be 
extrapolated to high energies, assuming power laws. 
The data are compatible with cut-offs in the spec- 
tra proportional to their nuclear charge Z . The all- 
particle spectrum and mean logarithmic mass thus 
obtained are compatible with air shower data. The 
combined energy spectra for groups of elements ob- 
tained from direct and indirect measurements ex- 
tend over seven orders of magnitude in energy. They 
are compatible with contemporary models of the ac- 
celeration of cosmic rays in supernova remnants and 
models of the diffusive propagation of the particles 
in the Galaxy [6] . 

Another important aspect is the anisotropy of 
the arrival direction of cosmic rays, giving hints to 
understand the propagation of the particles in the 
Galaxy. If the knee is caused by leakage of parti- 
cles from the Galaxy, one would expect less particles 



arriving from directions perpendicular to the galac- 
tic disk. Frequently, the Rayleigh formalism is ap- 
plied to quantify the observations [40] . Upper limits 
as obtained by KASCADE and preliminary results 
obtained by the KASCADE-Grande experiment are 
compiled together with results from many experi- 
ments in Fig. [5] [39] . An increase of the anisotropy 
amplitudes can be recognized. However, an inter- 
pretation has to take into account that the avail- 
able statistics decreases roughly as oc E~ 3 , there- 
fore, an increase of the observed anisotropy is ex- 
pected, which is of the order of the observed increase. 
The observed anisotropy amplitude is compatible 
with a diffusion model of the particle propagation 
[40)41(42] . 

5. TeV 7-ray astronomy 

Observation of TeV 7-rays with atmospheric 
Cerenkov telescopes, like the HESS [44], MAGIC 
[45] . or VERITAS (46] experiments, give the most 
direct hints towards the acceleration sites of cosmic 
rays. 

The energy density of galactic cosmic rays 
amounts to about pE ~ 1 eV/cm 3 . The power re- 
quired to sustain a constant cosmic-ray intensity 
can be estimated as L cr — p cr V/T esc w I0 41 erg/s, 
where r esc is the residence time of cosmic rays in a 
volume V (the Galaxy, or the galactic halo) . With a 
rate of about three supernovae per century in a typ- 
ical Galaxy the energy required could be provided 
by a small fraction (« 10%) of the kinetic energy 
released in supernova explosions. This has been re- 
alized already by Baade and Zwicky [47] . The actual 
mechanism of acceleration remained mysterious 
until Fermi [48] proposed a process that involved 
interaction of particles with large-scale magnetic 
fields in the Galaxy. Eventually, this lead to the 
currently accepted model of cosmic-ray acceleration 
by a first-order Fermi mechanism that operates in 
strong shock fronts which are powered by the explo- 
sions and propagate from the supernova remnant 
(SNR) into the interstellar medium [49150151152] , 

The theory of acceleration of (hadronic) cos- 
mic rays at shock fronts in supernova remnants 
is strongly supported by recent measurements of 
the HESS experiment [53)43] . observing TeV 7- 
rays from the shell type supernova remnant RX 
J1713. 7-3946, originating from a core collapse su- 
pernova of type Il/Ib. For the first time, a SNR 
could be spatially resolved in 7-rays and spectra 
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Fig. 7. The SNR RX J1713. 7-3946 as seen by the HESS experiment, spatially resolved in TeV 7-ray light. Left: 7-ray excess 
contours, superimposed are 14 fields for which spectra have been obtained independently. The corresponding integral flux is 
displayed as function of photon spectral index (right). For details see 1431 . 



have been derived directly at a potential cosmic-ray 
source as can be seen in Fig. [7J The figure shows on 
the left-hand side excess contours of TeV 7-rays. 
The SNR has been divided into 14 fields for which 
spectra have been obtained independently [43] . 
The corresponding integral flux values are plot- 
ted against the spectral index on the right-hand 
panel. The measurements yield a spectral index 
7 = —2.19 ± 0.09 ± 0.15 for the observed 7-ray flux 
averaged over the complete source. 

A model for the acceleration of hadronic particles 
in the SNR is supported by measurements in vari- 
ous energy ranges from radio wavelengths to TeV 7- 
rays [54] . An important feature of the model is that 
efficient production of nuclear cosmic rays leads to 
strong modifications of the shock with large down- 
stream magnetic fields (B ss 100 /J,G). Due to this 
field amplification the electrons are accelerated to 
lower maximum energy and for the same radio/x- 
ray flux less electrons are needed. Consequently, the 
inverse Compton and Bremsstrahlung fluxes are rel- 
atively low only. The model predicts that the spec- 
trum of decaying neutral pions, generated in inter- 
actions of accelerated hadrons with material in the 
vicinity of the source, is clearly dominant over elec- 
tromagnetic emission generated by inverse Compton 
effect and non-thermal Bremsstrahlung. The mea- 
surements are compatible with a nonlinear kinetic 
theory of cosmic-ray acceleration in supernova rem- 



nants and imply that this supernova remnant is an 
effective source of nuclear cosmic rays, where about 
10% of the mechanical explosion energy are con- 
verted into nuclear cosmic rays [54155] . 

Further quantitative evidence for the acceleration 
of hadrons in supernova remnants is provided by 
measurements of the HEGRA experiment [55] of 
TeV 7-rays from the SNR Cassiopeia A [57] and 
by measurements of the HESS experiment from the 
SNR "Vela Junior" (SNR RX J0852.0-4622) [55]. 

For completeness, it should also be mentioned 
that the detection of 7-rays with GeV energies con- 
tributed to the understanding of cosmic-ray propa- 
gation. The diffuse 7-ray background detected with 
the EGRET satellite experiment [59] exhibits a 
structure in the GeV region, which is interpreted as 
indication for the interaction of propagating cosmic 
rays with interstellar matter [60] , This is the most 
direct evidence for cosmic rays propagating in the 
galactic halo, well outside the galactic disc. 

6. Discussion and outlook 

In the last decade our understanding of the ori- 
gin of galactic cosmic rays has been significantly 
improved by multidisciplinary efforts, combining 
key observations of the direct and indirect mea- 
surements of cosmic rays as well as the detection 
of 7-rays. The observations by the HESS 7-ray 
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telescope give clear hints that hadronic particles 
are accelerated in SNR. The data are compatible 
with a model of first order Fermi acceleration at 
strong shock fronts. The particles propagate in a 
diffusive process through the Galaxy. Parameters 
of the propagation models have been constraint by 
direct measurements above the atmosphere. The 
KASCADE experiment has shown that the energy 
spectra of the light elements exhibit a cut-off struc- 
ture, while the spectra of heavier elemental groups 
follow power laws to higher energies. The observed 
spectra seem to be compatible with the assumption 
of power laws and a cut-off energy proportional to 
the nuclear charge. This implies that the knee in the 
all-particle energy spectrum is caused by a cut-off 
of the light elements. The shape of the all-particle 
spectrum at higher energies is then determined by 
the subsequent cut-offs of all elemental species in 
cosmic rays. Most likely, the astrophysical origin of 
the knee is a combination of the maximum energy 
reached in the acceleration process and leakage from 
the Galaxy during propagation [19]. More exotic 
ideas about the cause of the knee are most likely 
excluded [6l8] . 

In conclusion this gives a qualitative 'standard 
picture' of the origin of galactic cosmic rays. How- 
ever, several details remain unclear and a precise 
quantitative description of all aspects of the acceler- 
ation and propagation mechanisms is still missing. 
Among the open questions are: 
It is not clear how to precisely match the spectral in- 
dices observed at Earth to the spectra at the sources, 
being compatible with the TeV 7-ray observations 
and the modification of the spectral slope during 
propagation. 

A precise astrophysical interpretation of air shower 
data is at present limited by the understanding of 
hadronic interactions in the atmosphere, thus the 
exact shape of the energy spectra at their corre- 
sponding knees is unknown. 

Contemporary assumptions on the parameters of 
cosmic-ray propagation models yield anisotropies in 
the arrival directions not observed by air shower ex- 
periments [61119] . 

In the near future experiments like TRACER 
aim to reveal details of the cosmic-ray propagation 
for energies approaching the knee [62] . Experiments 
at the LHC probing the extreme forward direc- 
tion of phase space will improve the description 
of high-energy hadronic interactions [63] , Also the 
exploration of the end of the galactic component 
and the transition to extragalactic particles in the 



energy range from 10 17 — 10 18 eV will be of im- 
portance. Key experiments in this energy region 
are KASCADE-Grande (a 0.5 km 2 extension of the 
KASCADE experiment) [64], Ice Cube/ Ice Top (a 
1 km 2 air shower experiment and neutrino telescope 
at the South Pole) [65166] . and HEAT/ AMIGA (a 
25 km 2 extension of the Auger experiment to lower 
energies) [67] . 

Around 10 18 eV two features appear in the 
all-particle energy spectrum. The second knee at 
E2nd ~ 400 PcV« 92 x Ek, where the spectrum 
exhibits a steepening to 7 ~ —3.3, and the ankle at 
about 4 EeV, above this energy the spectrum seems 
to flatten again to 7 f» —2.7. The region around 
4 EeV is sometimes also called "the dip" in the 
spectrum. 

A possible cause for the second knee is the end of 
the galactic component, when all elements succes- 
sively have reached their cut-off energies, the latter 
being proportional to their nuclear charge |4|19| . 
If one assumes that ultra-heavy elements (heavier 
than iron) play an important role to understand 
the second knee, the factor of 92 between the ener- 
gies of the knee and the second knee can be easily 
understood as the nuclear charge of the heaviest 
elements in the periodic table. 
The dip is proposed to be caused by electron- 
positron pair production of cosmic rays on photons 
of the cosmic microwave background [68J . 

The investigation of the transition region and a 
precise measurement of the galactic all-particle spec- 
trum is also important for an estimate of the en- 
ergy content of extragalactic cosmic rays. For exam- 
ple, the extragalactic component needed according 
to the poly-gonato model to sustain the observed 
all-particle flux at highest energies has an energy 
density of p e ^ 9 = 3.7 ■ 10 -7 eV/cm 3 . 

A promising rediscovered technique for the ex- 
ploration of cosmic rays from the transition region 
to highest energies is the detection of radio sig- 
nals from air showers [69] . Most likely the emission 
mechanism is coherent synchrotron radiation of 
electrons with energies around the critical energy 
(85 MeV) deflected in the magnetic field of the 
Earth (geosynchrotron radiation) [70]. The LOPES 
experiment, registering showers in coincidence with 
the KASCADE-Grande experiment has demon- 
strated the feasibility of this approach [71] . Radio 
detection of air showers is also pursued in the CO- 
DALEMA experiment [72] and within the LOFAR 
radio telescope [73] . Also first radio pulses from air 
showers have been recorded with antennae set up at 
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the southern site of the Pierre Auger Observatory 
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